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The model of dense lattice polymers is studied as an example of non-unitary Conformal Field 
Theory (CFT) with c = — 2. "Antisymmetric" correlation functions of the model are proved to be 
given by the generalized Kirchhoff theorem. Continuous limit of the model is described by the free 
complex Grassmann field with null vacuum vector. The fundamental property of the Grassmann 
field and its twist field (both having non-positive conformal weights) is that they themselves suppress 
zero mode so that their correlation functions become non-trivial. The correlation functions of the 
fields with positive conformal weights are non-zero only in the presence of the Dirichlet operator 
that suppresses zero mode and imposes proper boundary conditions. 



05.50.+q, 11.25.H 

Introduction. — In spite of the remarkable success of 
unitary CFT's in predicting the critical properties of dif- 
ferent lattice spin models jjj, the non- unitary theories, 
although of no less importance for statistical physics, so 
far were not fully understood. It already becomes obvi- 
ous that some of the axioms of unitary CFT have to be 
sacrificed in this case |]. Still, it remains unclear where 
one has to modify the foundations and not to destroy the 
whole building of CFT. 

The general idea of the Letter is not to study the non- 
unitary CFT's on their own but, instead, to analyze one 
particular model of dense polymers on the lattice whose 
continuous limit corresponds to the non- unitary c = — 2 
CFT. We believe that at least some of the results ob- 
tained on this way should be universal and applicable to 
other non- unitary CFT's. 

The model of dense polymers actually has a long his- 
tory, dating back more then a century, when Kirchhoff 
proved a beautiful theorem that the number of one- 
components spanning trees (polymers) on the lattice of 
N sites is given by the principal minors of the N x N ma- 
trix of discrete Laplacian . Another fundamental re- 
sult was due to Fortuin and Kasteleyn ||[| who observed 
that the partition function Zn of the q-component Potts 
model can be represented as a dichromatic polynomial 
that continuously depend on q. Although the partition 
function of the model vanishes in the formal limit q — > 
owing to zero mode of the discrete Laplacian, its deriva- 
tive with respect to q does not and gives the partition 
function of one-component spanning trees. 

The purpose of the Letter is to show that: 

(i) The q — > limit of the Potts model can be car- 
ried on in two steps. The first, A — > 0, leads to 
the model of lattice polymers with arbitrary num- 
ber of components 7; the second, n — > 0, to their 
dense phase. Although the partition function of the 
model again vanishes in the limit, some "antisym- 
metric" 27-point correlation functions survive. 

(ii) These correlation functions are given exactly by the 
minors of rank (N — 7) of the Laplacian matrix. 



These can be rewritten in terms of integrals over 
anti-commuting variables and in continuous limit 
coincide with the correlation functions of the free 
complex Grassmann field. 

(iii) The vacuum vector of the field theory have to be de- 
fined as having zero norm. The fundamental prop- 
erty of the Grassmann field and its twist field (both 
are primary with non-positive conformal weights) 
is that their operator products define the Dirichlet 
operator that suppresses zero mode and imposes 
proper boundary conditions for the primary fields 
with positive conformal weights. This does not 
change other basic principles of CFT and leads to 
the logical and self-consistent theory. 

Dense Phase of Lattice Polymers. — Let lattice £ has N 
sites labeled 1,2, N. With each site i we associate a 
spin variable 0$ which can take q values, say 1, 2, q. 

Then the average of any operator A{o~) in the q- 
component Potts model we define as (without normal- 
ization factor!) 

(A(*)) = £ A(a) exp J fiJ^Sfavj) \ . (1) 

Here the c-summation is over all the spins o~±, er/v; the 
second summation is over all edges of the lattice. It has 
been shown that Zm can be expressed as a dichromatic 
polynomial To fix notations we briefly repeat the 

derivation of the result. Set v = exp(/3J) — 1, then the 
partition function can be rewritten as 

= <i> n I* ( 2 ) 

Let E be the number of edges of the lattice C. Then 
the summand in Eq.(Q) is a product of E factors. Each 
factor is the sum of two terms: 1 and v8((Ji,crj), so the 
product can be expanded as the sum of 2 E terms. 

Each of these 2 E terms can be associated with a bond- 
graph on the lattice C. To do this, note that the term is 
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the product of E factors, one for each edge. The factor 
for edge (ij) is either 1 or vS(cri, <Tj): if it is the former, 
leave the edge empty, if the later, place a bond on the 
edge. Do this for all edges (ij). We then have a one-to- 
one correspondence between bond-graphs on C and terms 
in the expansion of the product in Eq.(0). 

Consider a typical bond-graph Q, containing N sites, 
L bonds, 7 connected components and u> internal cycles. 
These are not independent, but must satisfy Euler's re- 
lation 



L + 7 = N + 1 



(3) 



Then the corresponding term in the expansion contains 
a factor v L , and the effect of delta functions is that all 
sites within a component must have the same spin a. 
Summing over all independent spins and over all bond- 
graphs Q that can be drawn on C we obtain |5||| 



Q 



(4) 



Note that here q need not be an integer. We can allow 
it to be any real number and, in particular, to consider 
formal limit q — > 0. Since we are going to deal with not 
only one- but arbitrary 7-component spanning trees, we 
have to treat the limit in a way different from || . 

At first we consider the limit X,q,v — > while n = q/\ 
and x = v/X remain finite. As a result we obtain the 
partition function of lattice polymers 

lim V k~<\"x l - ^ 



lim A 

A^O 



-A' 



z 



N 



+0 

g 



T 



K~>X L . 



(5) 



Zn = 



Here the last summation is over all bond-graphs T that 
has no internal cycles, i.e. lo = 0. Such graphs are usu- 
ally called spanning trees (polymers). The number of 
bonds L of the spanning tree is related to the number of 
its components 7 as L = N — 7. Hence, the partition 
function can be rewritten as 

= ^A/> 7 :r Ar -'', (6) 

7 Ti 7 

where symbol T 7 denotes the set of different 7-component 
spanning trees and A/" 7 is their total number. To simplify 
further notations we take x = 1 without loss of generality. 

The second limit n — > leads to the so-called dense 
phase of the polymer model. Since 7 > 1 the partition 
function (^|) obviously tends to zero in this limit. Never- 
theless, the correlation functions do not necessarily van- 
ish. Indeed, repeating all the steps leading to Eq. (||) one 
can calculate the following correlation functions 



lim (1) 

fc,A^0 

lim (Ski) 

s.A^O 



lim 

k,A^0 



Ski 

Sr,l 



= 0, 

= N(kl) = const, 

= N(kl)(pq) -N(kq)(pl), 



(7a) 
(7b) 

(7c) 



Here Ski — ( v /q)S{<Jk,&i)', N(kl) is the number of one- 
component spanning trees with both the sites k and I be- 
longing to the same component (this number, obviously, 
does not depend on the position of the sites); N(ki)(pq) is 
the number of two-component spanning trees with sites 
k,l belonging to one component and sites p, q to the 
other; etc. The antisymmetric combination of S's in each 
27-point correlation function is designed to guard against 
any contribution of spanning trees with the number of 
components less then 7 (otherwise this would be diver- 
gent). So, only 7-component spanning trees contribute 
to the 27-point correlation function in the limit k — > 0. 

The importance of these correlation functions is justi- 
fied by the following result. 

Generalized Kirchhoff Theorem. — Given a lattice C 
with N sites labeled 1,2, ...,7V, the N x N matrix of 
discrete Laplacian has the elements: =numbcr 
of edges incident to i, Ay = —number of edges with end 
points i and j. The minor A™' 1 ' of rank (N — 1) is ob- 
tained from the matrix A by deleting fc-th column and 
Z-th row; similarly, the minor 

is obtained by deleting columns k,p and rows l,q\ etc. 
Then 



det A 
detA«« 
det A^ 1 ^ 



0, 

: N(ki) 



(M)(pq) 



const, 



(kq)(pl), 



(8a) 
(8b) 
(8c) 



Here one immediately recognizes the "antisymmetric" 
correlation functions (0). The standard proof of the origi- 
nal version of the theorem (first two lines of the sequence) 
can be found in Ref. j|. Priezzhev [Q proposed an al- 
ternative proof of the original version in the spirit of the 
combinatorial solution of Ising model. His method is sim- 
pler and can also be generalized to prove all other lines 
of the sequence (||). 

Free Complex Grassmann Field. — Using the matrix 
representation we can reinterpret the partition function 
of lattice polymers as being the partition function of some 
artificial statistical system. To this end we define at each 
site i of the lattice C the pair of anti-commuting variables 
9i and 9* (its complex conjugate). Then, using Berezin's 
definition of the integral over anti-commuting variables 
|Q we can rewrite the determinant of the matrix A as 

det A = f d9{ ...d9 N exp ^ 9* & lJ 9 3 (9) 

= f d9{ ...d9 N exp ( 9 *i Oj ) ■ 

In continuous limit this partition function defines field 
theory with the action 



S[9] = -!- / d^trt) d 2 r. 



(10) 
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The average of any operator A[9] we define as 



(A[6]) 



~d9*d9] A[0] exp-S[0] 



(11) 



Then, although the average of the identity operator is 
equal to zero due to the presence of zero mode, all other 
correlation functions of the field 9 are non-trivial and can 
be normalized so that 



(1 



= 0, 
= 1, 

= ln(% 1 4 2 ) : 
m(r?i) 



ln(% 2 l) 



(12a) 
(12b) 
(12c) 

(12d) 



The field 6 is scalar and its correlation functions depend 
only on the projectively invariant cross-ratios 



nil 



\r14r23) 



(13) 



Here 61 = 9(ri); n.2 = \v\ — v%\. These correlation func- 
tions are nothing but asymptotics of the "antisymmetric" 
correlation functions ([[]) in the continuous limit. 

The surprising thing is that the Grassmann field itself 
suppresses zero mode of the Laplacian operator. In spite 
of this unusual property it still can be considered as a 
primary conformal field with the weight hg = 0. 

The stress-energy tensor, 



T(z) 



-.89*89: 



lim <^ 89*{z)89(w) 



satisfies standard operator product expansion, 
-1 2T(w) dT{w) 



T(z)T(w) = 



w 



(14) 



(15) 



One can assure himself by direct calculation that the 
stress-energy tensor is indeed the generator of confor- 
mal transformations in the sense that for any correlation 
function (X) = (0* . . . 9 2 n) from the sequence ( ]l2| ) its 
transformation law is given by 

S e {X) =<£ dz e{z) (T(z)X) + f dz e(z) (f(z)X). (16) 
Jc Jc 

The correlation functions of the field 9 satisfy the third- 
order differential equation coming from the condition of 
degeneration of the operator (1, 3) with the weight hi 3 = 
on the third level. This equation actually becomes of 
the second order for the field 89 and, in its turn, coincides 
with the condition of degeneration of the operator (2, 1) 
with the weight h2,i = 1 on the second level. 



The twist field cr(z, z) can be defined with the use of 
the standard operator product expansion 



89(z)a(w, w) 



t{w, w) 



w 



(17) 



Alternatively, on the lattice it can be defined by means of 
the construction similar to that for the disorder operator 
in Ising model || . 

Conformal properties of the twist field a are similar 
to those of the Grassmann field 9. Namely, its correla- 
tion functions are non-trivial even in the presence of zero 
mode. Its correlation functions can be found from the 
condition of degeneration of the operator (1, 2) with the 
weight hi 2 = —1/8 on the second level || 



(aia2&3&4:) 



'7*12, 



7TV r 12 r 34\/| r ?(l 

{F( V )F(l - fj) - 



F(v)F(l 



(18a) 



77)}, (18b) 



where F(ri) = 2^1(5, \\ and V = (^13^24) / '(z 12 z 34 ) . 
Mixed four-point correlation function of the fields 9 and 
a can also be found using standard techniques of CFT 



{0*i6 2 a 3 a 4 ) = 2^{H( V )+H(fj)} 



(19) 



Here H{rf) = In (y 7 ^ + y 7 ?; — l) . 

This means that both the Grassmann field 9 and its 
twist field a can be considered as primary conformal fields 
with the weights hg = and h„ = —1/8 provided that 
the vacuum state has been defined as having zero norm. 
These fields are unique in having both the property and 
non-positive conformal weights. 

Dirichlet Operator and Green Function. — There is a 
simpl e re lation between the four-point correlation func- 
tion ( 12c ) and the Green function of the Laplacian oper- 
ator. The most straightforward way to understand this 
is follows. Let us consider a conducting plane with a cur- 
rent 1=1 entering the plane at a point ri and leaving it 
at a point r 2 . Then the voltage difference between sites 
r 3 and on the plane is given by the four-point function 

The Green function of the Laplacian operator with the 
Dirichlet boundary conditions at the point ro can be de- 
fined quite similarly. Consider the same conducting plane 
earthed at the point ro. This means that the voltage at 
this point is always maintained to be equal to zero. If a 
current 1=1 enters the plane at a point ri (and leaves 
it at the earthed point ro) then the voltage at a site r 2 
is given by the Green function Go( r i> r 2 ). 

The operator T>q that corresponds to the earthed point 
ro can, obviously, be considered as the product of the 
field #0 with its complex conjugate #0 a t t ne same point. 
We will call it the Dirichlet operator since it imposes the 
Dirichlet boundary conditions on the Grassmann field. 
With the help of this operator the Green function can be 
represented as 
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G (ri,r 2 ) = {V Q 6ie 2 ). 



(20) 



The Dirichlet operator can formally be denned through 
the following operator products 



V = lim{6>n6>i} = lim 



Vnii 



(21) 



Its two-point correlation function can be found from the 
four-point functions (12c,18t .1^). However, one has to be 
careful merging different points of the four-point func- 
tions because they diverge logarithmically in the limit. 
These divergences have absolutely the same nature as 
those present in the Green function in the thermody- 
namic or continuous limit To treat them carefully 
let us first consider the correlation functions on the lat- 
tice with spacing a. This scale dictates minimal possible 
distance between different merging points. After arbi- 
trary conformal transformation the lattice is no longer 
uniform and the area of any given fundamental square of 
the lattice acquires an additional factor proportional to 
the metric on the plane: ds 2 = <?(r)dr 2 . Finally, we have 



(the factor 



Pi 



is absorbed into the metric!) 



1, 



(T>iT>2 



In 



(n 2 ) 4 
9i 92 



(22) 



Note, that the Dirichlet operator is scalar and its correla- 
tion functions are projectively invariant. This is natural 
since it has been defined as the product of scalar fields 
9 and 9*. This also suits its interpretation as being the 
operator that determines boundary conditions. 

The operator product of the Dirichlet operator with it- 
self and with the operators 9 and a cannot be completely 
determined within the CFT. Indeed, merging different 



points of the correlation functions (12) we can only say 
that 



lim {0„2M 

l-i-O 


= k 9 , 


(23a) 


lim {ct X>i} 


= k cr , 


(23b) 


lim {2?oZ>i} 

l->0 


= kV , 


(23c) 



where k is some constant. From the point of view of the 
"electrical" interpretation of the operator given above the 
most natural choice would be k = 1. 

As an example of the field with positive conformal 
weight let us consider correlation functions of the local 
energy operator 



e = .8^9*8^9: = Km{dvfodPOi - 47r<3(r i)}. 



(24) 



This is primary with conformal weight h e — 1. Its cor- 
relation functions can be found from Eqs. ( |l2| ) and are 
all trivial, (e\ . . .en) — 0, unless we insert the Dirichlet 
operator 



(X> £l£2) 
(X> £l£2£3£4) 



64 
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G4 



(ri 2 r 3i ) 4 (ri 3 r 2 4) 4 (ri 4 r 2 3) 4 



(25a) 
(25b) 



This property is common to all primary operators with 
positive conformal weights. 

Let us summarize the results of the Letter. It has been 
shown that the model of the free complex Grassmann 
field properly describes the continuous limit of the lattice 
model of dense polymers only provided its vacuum vector 
has been defined as having zero norm, (0|0) = 0. Never- 
theless, it is this vacuum that has to be considered when 
one studies the correlation functions of the primary fields 
with non-positive conformal weights (9 and a). These 
correlation functions, ((l^), imply the mode expansion 

6{z, z)= xa + 20o ln\z\ - (—z~ n + — *" n ) , (26) 

n^o ^ n n J 

with the commutation relations 

{Xo, Xo} = 3. {6*, 9 m } = {9* n , 9 m } = n S n+m , (27) 

where (3) = 1 and (Q 2 ) = 0. Operator 9 is nothing but 
the coordinate-independent part of the Dirichlet opera- 
tor. It defines yet another null vector, |*) = 3|0). This 
can be normalized so that (*|0) = 1. Together these 
two vectors, |0) and |*), define physical vacuum state for 
those primary fields that have positive conformal weights. 

We conclude that the theory is non-trivial only due 
to the presence of two different null vectors that are not 
orthogonal to each other. This could be a general feature 
of other non- unitary CFT's. 

I would like to thank V.B. Priezzhev for many stimu- 
lating discussions. 

Remark. — After this Letter was written, V. Gurarie 
and M. Flohr informed me that the first two lines of the 



sequence (12i,b) have already appeared in their paper |9[| 
where conformal properties of the Quantum Hall State 
have been analyzed. They also considered the product of 
the fields 9 and 9* and noticed that the resulting operator 
suppressed zero mode. However, there is an important 
difference between the results of Ref. (see also ear- 
lier references therein) and this Letter. Namely, in this 
Letter it has been shown that for the CFT to be consis- 
tent with the lattice model of dense polymers the field 8 
have to be considered as a primary field with the weight 
he = 0. Hence, its product with its complex conjugate 
(the Dirichlet operator) has to be invariant under the pro- 
jective transformations SL(2,C). This is important since 
it is this operator that determines vacuum state for the 
fields with positive conformal weights. On the contrary, 
the authors of Ref. || did not consider the field 9 as 
conformal and omitted the metric factors in Eq. (|2^). 
The "logarithmic" operator they obtained in this way 
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was not scalar but possessed more complicated transfor- 
mation properties. 
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